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Propagation of three dimensional magnetosonic waves is considered for a homogeneous shear 
flow of an incompressible fluid. The analytical solutions for all magnetohydrodynamic variables 
are presented by confluent Heun functions. The problem is reduced to finding a solution of an 
effective Schrodinger equation. The amplification of slow magnetosonic waves is analyzed in great 
details. A simple formula for the amplification coefficient is derived. The velocity shear primarily 
affects the incompressible limit of slow magnetosonic waves. The amplification is very strong for 
slow magnetosonic waves in the long-wavelength limit. It is demonstrated that the amplification of 
those waves leads to amplification of turbulence. The phenomenology of Shakura-Sunyaev for the 
friction in accretion disks is derived in the framework of the Kolmogorov turbulence. The presented 
findings may be the key to explaining the anomalous plasma heating responsible for the luminosity 
■ of quasars. It is suggested that wave amplification is the keystone of the self-sustained turbulence 

(N ■ in accretion disks. 
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• , I. INTRODUCTION 

I ■ The occurrence of intense dissipation in accretion flows is among the long-standing unsolved problems in 

O I astrophysics^^ Revealing how the magnetized turbulence creates shear stress tensor is of primary importance to 

■^-^ I understand the heating mechanism and the transport of angular momentum in accretion disks. The transport of 

' angular momentum at greatly enhanced rates is important for the main problem of cosmogony, that is understand- 
ing the dynamics of creation of compact astrophysical objectsi^ Without a theory explaining the enhanced energy 
dissipation in accretion flows of turbulent magnetized plasma we would have no clear picture of how our solar system 

I has been created, why the angular momentum of the Sun is only 2% of the angular momentum of solar system, while 

' carrying 99% of the solar system's mass, why quasars are the most luminous sources in the universe. The importance 

QQ I of friction forces and convection as well the problem of angular momentum redistribution for the first time was em- 

. phasized by von Weitzacker^^ a very detailed bibliography on the physics of disks is provided in the monograph by 

' Morozov and Khopcrskovi^ Gravitational forces, angular momentum conservation, and dissipation processes become 

^Kf^ , the main ingredients of the standard model of accretion disks by Shakura and Sunyae"\*i and Lyndcn-BcU and Pringlci^ 

' Lynden-Bell^ suggested that quasars are accretion disks and Shakura and Sunyaev introduced alpha phenomenology 

Q^ i for the stress tensor 

o 



(^RV^OLP^ P'^Pf^s^ (1) 



■ where a is a dimensionless parameter, p is the pressure p is the mass density, Cg is the sound speed, and the indices 
5—1 ' of the tensor come from the cylindrical coordinate system {R, (p, z) related to the disk rotating around the z-axis. 
As accretion disks can have completely different scales for protostellar disks, mass transfer disks, and disks in active 
galactic nuclei (AGN) it is unlikely that Coriolis force is the main cause for dissipation, while the bending of the 
trajectories is a critical ingredient. We suppose that the shear dissipation in magnetized turbulent plasma is a robust 
and very general phenomenon which can be analyzed as a local heating for approximately homogeneous magnetic field 
and gradient of the velocity. For these reasons in the next section we will use local Cartesian coordinates choosing 
the X-axis in radial direction e^. = R, and the j/-axis along the circulation of the almost Kcplcrian motion of accreting 
plasma, ej^ = ip. The axis of the disk is along = z. 

There is almost a consensus that the magnetic field is essential and should be introduced from the very beginning in 
the magnetohydrodynamic (MHD) analysis. The likely importance of MHD waves on the cosmogony of solar system 
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has been pointed out by Alfven: "At last some remarks are made about the transfer of momentum from the Sun to 
the planets, which is fundamental to the theory. The importance of the magnetohydrodynamic waves in this respect 
is pointed out.—" As the magnetic field lines are frozen in the highly conducting hot plasma, we suppose that they lie 
in the plane of the disk B « —Be^p. In our local analysis we can suppose that mass density is also constant p « const. 
Moreover, we assume that the magnetic field is sufficiently small and the Alfven speed is much smaller than the sound 
speed, Va ^ Cs. As usual, the magnetic pressure is defined as 

\pVl=pB^^. (2) 

Under these conditions, we can consider the plasma as an incompressible fluid. In other words we assume a plasma 
for which the magnetic pressure ps is much smaller than the gaseous one p 

/3 = — » 1. (3) 

PB 

According to the theory of turbulent dynamo the density of magnetic energy ps is comparable with the density of 
kinetic energy of the turbulent motion at a basic scale which is smaller than the pressure, i.e., we suppose that the 
turbulence creating the magnetic field is subsonic. 

Our goal is to demonstrate that small magnetic fields can catalyze a disk to ignite as a star and to help the kinetic 
energy of the shear flow to be converted into heat in spite of the very low molecular plasma viscosity rj. 

In a weak magnetic field the incompressible transverse MHD waves, the slow magnetosonic waves (SMWs) and 
Alfven waves (AWs) , are qualitatively new features of the weakly magnetized plasma. Under these conditions every 
mechanism of giant wave amplification inevitably acts as a dissipation mechanism for heat production which creates 
stress tensor and a transport of angular momentum in accretion disks. 

The present research is triggered by the amplification of SMWs observed by numerical analysis of two-dimensional 
(2D) MHD wavesjii Even these first investigations demonstrated that shear flow leads to exchange of energy and wave 
amplification, mutual transformation between different wave modes, gave the perspectives to explain self-heating 
and other processes now known as nonmodal. As we have standing wave amplification often is used the notion 
overre flection] a running wave with ky is amplified together with the reflected wave with wave- vector —ky. These 
conclusions was confirmed extended and popularized by the numerical investigationsii^ However without a realistic 
analytical three-dimensional solution the investigations of shear flows is in state of infancy. The purpose of our work is 
to analytically solve the simplest case of three-dimensional (3D) waves in homogeneous shear flow and magnetic field, 
to analytically calculate the phase averaged amplification coefficient A of the waves and to incorporate this waves's 
amplification into some standard model for turbulence. Or, in short, our final aim is to use an exact solution for 
linearized waves in order to build up a coherent scenario for accretion disk theory and luminosity of quasars. Before 
continuing we wish to refer to the concluding remarks from the review by Balbus and Hawley>^ Again and again we 
are ignorant. The good news is that, for first time, it appears that we know in which directions we should be looking 
to begin to find answers to questions like these. 



II. MHD MODEL. ANALYTICAL SOLUTION TO LINEARIZED WAVE EQUATIONS 

Our model problem is to investigate MHD waves in a homogeneous shear flow and a homogeneous magnetic field 

Vshciir ^ Axey, Bo = -Soe.y, e, = d^r, eg = Bq/Bo = -e^, (4) 
= V X V.hear = Ae^, i = 1, 2, 3, eo = n/n = (5) 

in an ideal and incompressible plasma. For the linearized MHD equations we suppose a time dependent wave-vector 

k(t) = —AtkyBx + kyBy + k z , k ' T = (y — Atx)ky + k^z (6) 

in the plane waves solution. For the velocity V and the magnetic field B we consider small enough time-dependent 
amplitudes v and b 

V = Vshcar + Vwavo , B = Bq "f B^avo , (7) 

Vwavc = -il^Av(t)e* ^ Bwavc = B^\y(t)(P^ \ 

The magnetic field B and velocity V are, of course, real vector fields as functions of t and the position vector 
r ~ (a;, y, z\ but it is mathematically convenient to use complex variables and calculate the real parts of the derived 
solutions. 
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Let us summarize the supposed approximations: incompressible V • V = 0, and inviscid 77 = 0, fluid with constant 
density p = const, and negligible Ohmic resistivity p = 0. For a sufficiently slow nonrelativistic motion the MHD 
equations take the formic 

pDtV = -Vp - B X (V X B)/^o, DfB = B • VV, (8) 

where Di = 9f + V • V is the substantial Lagrange derivative; V • B = 0. 

Let us introduce the characteristic length-scale A and velocity Va of the problem 

A=-?^<i?, Va^^^, 1^a=|Va|«Cs, (9) 



dimensionless wave- vector K = Ak, and dimensionless time 

KM) kJt) , , 

T = At= '±1 = — ^A_L_ 10) 

y y 

The initial time is chosen so that k^it = 0) = 0. Wc suppose that the characteristic length A is much smaller than 
the disk thickness ddisk- 

As the shear parameter describes the radial gradient of the orbiting velocity 

^ = - -2^Kcplor (.ttj 

at the corresponding radius R, it is of order of the angular frequency ti^Kopier of the Keplerian motion. We suppose also 
that for weak magnetic fields the Alfven speed is much smaller than the velocity of Keplerian orbiting, Va ^ Vxcpicr = 
RujKcpioY- In this case the wavelength of the slow magnetosonic waves is much smaller than the radius of the orbit 

R\Jk'^ + ^ 1 and the disk thickness Ifc^j ddisk ^ 1 and we approximately consider the shear flow as homogeneous. 

In other words, the goal of the present research is to investigate the statistical properties of a shear flow determined by 
short wavelength waves inaccessible by direct astronomical observations and usual computer simulations of accretion 
disk magnetohydrodynamics. 

The Coriolis force has a negligible influence on the propagation of magnctohydrodynamic waves and their amplifica- 
tion analyzed in the present work for the case of dominating magnetic field. Consequently, in our further consideration 
we will neglect that force. In such a way we have to investigate the influence of a homogeneous shear flow on small 
amplitude MHD waves in an incompressible ideal fluid. The homogeneous magnetic field we suppose to be parallel 
to the shear fiow. 

With so introduced variables the substitution of wave components from Eq. ([7]) in MHD equations Eq. ([8]) after 
some algebra yieldsii 

AtVx = 2 ^ "^^ ^ Kybx, <^Tbx = KyVx, (12) 

.KyK_Z 

K-v = = K-b, K = -KyTex+Kyey + K;,e^. (14) 



drV^ = 2—^p^Vx ~ Kyb^, d^6z = KyV^, (13) 



An unabridged derivation of this set of equations is given in Ref. [l5|. Later on, the substitution of drV^ with <£^hx 
in Eq. and analogous substitution of drVz in Eq. (|13p leads to 



1+e 



dp, + Q% = 2K,^^, (16) 



where 



Ky Kx 

^^T-f = -^^=, (17) 



Q = JK^ + Kl = const. (18) 
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For a two-dimensional motion {K^ = 0) in the plane of the disk £, = t. The negative friction for ^ < in the oscillator 
equation Eq. p5p comes from the time-dependence of the x-component of the wave- vector 

K.,{T)=-KyT = -Qt (19) 

The other components are constant Ky = const, and = const. In such a way the change of the sign in the effective 
friction force is a property of the shear flow related to the change of the sign of the corresponding component of the 
wave- vector. 

For the Alfven wave amplitude in Eq. (|16p we have the equation of a harmonic oscillator with an external force. 
This equation has the Green function 



(d| + Q'-) Git a = m - e') = diOi^ - e'), (20) 

Q 



a^^^^^^^^mti^oi^.a (21) 



and the general solution to Eq. (|16p reads 



(d| + g') xiO = 0, xiO = Qxg + C^Xn, 

Xg(0 = cos(QC), ;^,(0 = 5^^, (22) 

where Cg and Cu are the amplitudes the even and the odd eigen ^-oscillations related to the pure Alfven waves (AWs) 
for the case of = 0. AWs are not influenced by the shear flow. For the solution of the Cauchy problem with fixed 
bzi^o) and d^bz{£,o) = Qvz{0) we have 

-^l' MQi^ " 0]Y^dt + b.Mo) cos[Q(C - Co)] + d,54^o)^^^^i^|^. (23) 

For the special case of = we have a complete separation of variables in two independent sets of equations. One 
of them is for Alfven waves with oscillating and Vz perpendicular to the wave-vector and magnetic field, and the 
second set is for the independent slow magnetosonic waves (SMWs) for which the oscillations of the magnetic field 
and velocity are in the plane of the external magnetic field and the wave-vector. Still we can use the terminology 
AWs and SMWs for waves having nonzero Kz- In this case the amplitude of the SMW bx plays the role of external 
driving force for the amplitude of AW in Eq. . We will see later that amplification is related to SMWs, and AWs 
with bz components are only influenced by SMWs at nonzero Kz- Due to a resonance this influence, however, could 
be significant. It is notable that AW amplitude bz has provides no feedback on the SMW amplitude 6^. 

MHD equations give the possibility to express the velocity by the magnetic field. For the velocity from Eq. 
and Eq. ([13]) we have 

The incompressibility conditions Eq. p4p yield explicit expressions for the oscillations parallel to the external magnetic 
field and shear flow 

K^b^ + Kzbz -KyTbx+Kzbz Q K^v^ + KzVz 

by = = ^—^bx-—bz, vy = ^^v^- Vz. (25) 

y y y y y y y 

In such a way the main detail in solving the set of Eq. p2p and Eq. p3|) is to derive an analytic formula for b^. 
The first derivative A^b^ term in Eq. (jl5[) disappears by the substitution 

and we arrive at an effective Schrodinger equationi^ which describes the SMW 



dii^ + 



1 



il) = Q. (27) 
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This "Schrodingerization" gives the possibihty to use the quantum mechanical analogiesi^ and also the whole system 
of notions of quantum scattering theory. Constant Wronskians give also a technical convenience. As we will see later 
the amplification of SMW in 2D case {k^ = 0) is represented by the asymptotics phases of the wavefunctions Eq. ([88|l: 
a simple illustration of Heisenberg 1938 idea that S-matrix contains whole experimentally accessible information. 

As the effective potential 2mU /h^ = 1/(1 + ^^)^ is an even function of ^, the general solution to the Schrodinger 
equation Eq. ()27|1 is a linear combination of the even and odd solutions 

^(0 = Cg^g(C) + a^u(0, (28) 

which obey the boundary conditions 

Vg(0) = 1, djVgW = 0, (29) 

Vu(0) = 0, A^M^) = 1. (30) 

The analytical solutions are represented by the confluent Heun function 

V'g=v/rTe^HcunC(0,-i 0,-^,i±^,-a, (31) 

^u=eyT+C^HeunC(0, +1,0, -^l±^^^e)- (32) 

In the computer algebra system Maple the confluent Heun function reference is HcunC(a, /3, 7, (5, 77, z). This special 
function obeys the equation 

z{z - l)y" + [Az^ + Bz + C] y' + [Dz + E] y = 0, (33) 

where 

A = a, B = 2 + /3 + 7-a, C = -l-/3, (34) 

D = i[(2 + 7 + /3)a + 2(5], E = i[-a(l + /3) + (1 + 7)/? + 7 + 277)]. (35) 



For the series expansion 



we arrive at the recursion 



where 



y = ^c„ = ^a„z", c„ = a„z" (36) 



n=0 



Fa„_i + Ga„ + Ha„+i = 0, (37) 



F = 5+i(2n + /37), (38) 
G = n2 + (l + 7 + /3 - a)n + i [7 + 2r/ - a + /3(1 - a + 7)] 

H= -(n+l)(n+l + /?). (39) 
Supposing a_i = and oq = 1 we use the recursion 

a„+i = -(Fa„_i + Ga„)/H (40) 

which, for example, gives 

ai = i[/3(l + 7-a)+7 + 2r/-a]/(l + /?). (41) 

As the effective Schrodinger equation Eq. (j27p has a solution for arbitrary ^ € (—00, +00), the formal scries for Heun 
function have convergent Fade approximants. 
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Those Pade approximants can be calculated by the well-known e-algorithm. First we calculate the series of the 
partial sums in zeroth (0) approximation 



= Co + ci + 



c(0) _ c(0) , 

— "I" Ci 

and for the first N terms we get 

Sq = Co, — Co + ci, . . . , 

Our problem is to calculate the limit of the sequence 

S = lim Sn- 

n — >oo 

For this calculation we generate the auxiliary sequence 

Hi"' = l/c.+r 

recalling that 1/0 = 0, i.e., using pseudoinverse numbers if we have to divide by zero: 



+ CAT. 



hI,"' = l/ci. 



h["^ = l/c„ 



H 



(0) 
N-l 



I/CN- 



The epsilon-algorithm is the calculation of the recursion for the series 



s: 



(fc) 



4+1 



H, 



(fe) 



(fe-i) 
i+i 



5, 



(fe-i) 



(42) 
(43) 
(44) 

(45) 

(46) 

(47) 
(48) 



I+l 



gives the best Fade approximant 5*1^^ for the searched 



for all indices for which these relations make sense. 

The maximal in modulus auxiliary element i/max = H 

limes and the accuracy is of the order of l/i?max- 

For each accuracy of the final result e we can calculate S and H sequences with an accuracy S to assure that 
l/i/max < £• In such a way we obtain a method for calculating the confluent Heuu function and the solution to the 
MHD equation in power series of time. If from physical arguments we know that a solution exists, the divergent series 
can be summed. Having a method for calculating the confluent Heun function, we can calculate both the even and 
odd solutions to Eq. The accuracy in calculating Heun functions is controlled by the Wronskian 



w^(V'g,V'u) 



MO MO 
diMO d^MO 



= 1. 



(49) 



The constants from the general solution are also given by the Wronskians 

Cg = W^(^, Vu) = Mo)diAMo) - M^o)dM^o), Cu = W^(^g,^) = M^o)dM^o) - Mo)diM^o). (50) 

Those formulae generally apply for a Cauchi problem where the initial conditions are imposed on the function being 
sought, i.e., on V'(Co) and its derivative dj'0g(?o)- 

All MHD variables can be expressed by the solution Eq. (pS)) to the effective Schrodinger equation. Then bx{^) is 
given by Eq. and bz{^) by Eq. (f^^ . Equation (|25p gives by. And finally we know the time t = Q£^/Ky dependent 
amplitude of the magnetic field b(^ = Kyr/Q). Analogously the velocity can be expressed by Eq. p2|) 



V(e) = dQ5b(0, dg^ 

i.e., plasma's displacement is parameterized by the magnetic field 



QdC' 



b(0 = Q j ^rdi^VAky j ^^At 



(51) 



(52) 



recalling that magnetic field lines are frozen into the highly conductive fluid. 
The physical time t is related to the dimensionless one r 



t = 



A AK., 



(53) 



y 
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The time-dependence of the wave- vector Eq. ([6]) allows us to express the magnetic field B(t, r) and velocity of the 
fluid V(t, r) Eq. ([7]). The general solution depends on 4 arbitrary constants Cg, Cu, Cg, Cu- 
After some algebra we can express also the wave component of the pressure^^ 

Pit, r) = pV^Pit) c'-'W- ^ Pit) = -hyit) + 2^1^. (54) 

In such a way we derive a general analytical solution for linearized MHD waves in a shear flow of a magnetized fluid. 
In the next section, we will analyze the behavior of this exact solution when the waves amplification is significant. 

The essential part of our analytical solution for three-dimensional SMWs in a shear flow is presented by the x- 
component of the magnetic field 

-B, =BoRe{exp[i(k(i) •r + y.o)]} (55) 
X |cgHeunC(0,-i 0,-^,i±^,-e^) 

+Cu eHeunC(0, +1,0, -e^)| , 

= Bo cos [k(i) •r + (po]{Cg6^,g(0 + Cu6^,u(C)}, 

where 

g^^r= Ait-to), (56) 

Q' = K\kl + kll = (57) 

Here lpq and to a-i'c arbitrary constants parameterizing the initial conditions. All other MHD variables: velocity, 
pressure and displacement of the plasma by the wave can be expressed by this solution. The dimensionless times have 
simple geometrical interpretation in the k-space 

k^it) k.,{t) 
k — ' / (5«) 

For the velocity of AWs a substitution of 4' from Eq. ((26l) into Eq. ([22|l gives 

v.iO = dQ,h = -Cg smiQO + a^^^ + ^ cos[Q(C - C')] (59) 



If we substitute here ip from Eqs. (|28p. (|3ip and (|32|) . we arrive at the final analytical solution which contains Hcun 
functions. 



III. AMPLIFICATION OF THE MHD WAVES 



The odd and even solutions have asymptotics at ^ ^ oo 

V'g « IJg cos(QC -I- Sg), ?/'u « i'u cos(Q^ + Su), (60) 

where the asymptotic phase shifts 6giQ^), S^iQ^) depend on the effective energy. For a sufhciently large wave- vectors 
we have the asymptotic 

<5g(Q2 » 1) ^ 0, 6^iQ'- » 1) « -|. (61) 

As we will see later, the averaged amplification coefficient of the energy of MHD waves AiSg, S^) depends only on 
the asymptotic phases of the solutions but not on the amplitudes Dg and I?u- 
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A. Auxiliary Quantum Mechanical Problem 



Temporarily introducing imaginary exponents exp(ik • r) instead of sin(k • r) and cos(k • r) gives significant sim- 
plification of the analytical calculations related to the physics of the waves. In this subsection we will consider 7/; 
in Eq. ([77]) as a complex function in order to make easier any further analysis of the MHD amplification coefficient. 
In order to analyze the effective MHD equation Eq. ((27)) wc will solve the quantum-mechanical counterpart of our 
MHD problem, that is a tunneling through a barrier 2mU/fi^ = 1/(1 + C^)^, supposing that ip is a complex function. 
Consider an incident wave with a unit amplitude, a reflected wave with amplitude R and a transmitted wave with 
amplitude T 

iPiC -00) « exp(iQO i?exp(-iQO, (62) 
V'(C ^ +00) « Tcxp{+iQO. (63) 

Using the asymptotics of the cigcn-functions 

^ / cos(Qe - Sg) for e ^ -00, , 

"^^^ \Dg cos(Qe + Sg) for ^ ^ +00, ^^^> 

and 

_ / -Dn cos(gC - 5u) for ^ -00, , . 

V'^- \ cos(Qe + for e ^ +00 

as well as the general condition 

m = C^^''UO + C\^Hn{S.), (66) 

we compare the coefficients in front of exp(i(5^) and exp(— iQ^) for ^ —00 and ^ — > 00. The solution to a simple 
matrix problem yields 

C(^)i?g = exp(i5g), Ci'^^D^ = - exp(i,5u) (67) 
Then for the tunneling amplitude we get 

T ^ |r|ei('5"+'5s-'^/2) (68) 
and finally for the tunneling coefficient we obtain 

V = |T|2 = s^g, sug = sin((5u - <5g). (69) 

The convenience of the tunneling coefficient is that it varies in the range < P < 1. In the next two subsections we 
will present the SMW amplification coefficient as a function of the tunneling coefficient /C = 2/2? — 1 on the analogy 
of Heisenberg's ideas in quantum mechanics where the statistical properties of the scattering problem depend only on 
the phases of the S'-matrix. 

The phases dg, 6^ & (— tt, tt) and amplitudes Dg and can be determined continuing the exact wave functions 
Eqs. with WKB asymptotics Eq. Thus we obtain 

i; = D cos(Qe + (5 = ~QC - arctan , 

D= (70) 

cos(Qe + ^)' ^""^ 

at some sufficiently large ^ 3> 1 + 2tt/Q. Here int(. . . ) stands for the integer part of a real number. When programming 
we have to use the two-argument arctan function 

TT 

arctan(?/,x) = arctan(2//2;) + — 9{~x) sgn{y) £ (— tt, tt). (71) 

The accuracy of this continuation is controlled by the Wronskian from the asymptotic wave functions 

Wiijg, = QDgD^ sin(5g - Su) = 1. (72) 
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B. MHD and Real i/' 

Imagine that in an ideal plasma we have at t ^ — oo some plane MHD wave ~ our task is to calculate how many 
times the energy density increases at t — > oo, and to average this amplification over all the initial phases of that wave. 
As there is no amplification for the component according to Eq. we will concentrate our attention on the 
component. The amplification comes from the negative "friction" term oc ^/ (1 +^^) in Eq. (fT5l) . The influence of this 
friction is transmitted to the effective potential barrier oc 1/(1 + ^^). For < 1 we have an analog of the quantum 
mechanical tunneling. 

In the current MHD problem ip is a real variable with asymptotics 

~ / cos(QC - (f^i) for ( ^ -oo, , . 

'''^ \Dicos{Q^ + (f>i) fore^+oo, ^'"^^ 

i.e., we have an incident wave with a unit amplitude and an initial phase Df is the amplitude and is the phase 
of the amplified wave. 

Again we present the ip function as linear combination of even and odd solutions 

^'g(C)+C(^Vu(C). (74) 

Here we substitute the asymptotic formulas Eq. ([64|) and Eq. (|65|l . and the comparison of the coefficient with Eq. ([73|) 
at ^ ^ — oo gives 

C(^)D^=^^, Ci^)D^ = ^, (75) 

Sgu Sgu 

where 

Siu = sin((/)i - (5u), Sig = sin(0i - (5g). (76) 
The comparison of the coefficients at ^ ^ oo gives for the phase and the amplification of the signal 

SigSu ~t~ SiuSg 
SigCu -t- SiuCg 



Df = F(0i) = arctan ^'g"" ' , (77) 



where 



A{<P,)^D^ = ^, (78) 

^f = (SigSu + SiuSg)^ + (SigCu + SiuCg)^, (79) 

Sg = sin((5g), Su = sin((Su), (80) 
Cg = COs(^g), Cu = cos((5u). (81) 

The reversibility of the dissipation-free motion leads us to 0i = F{4i{)^ i.e., function F coincides with its inverse 
function F{F{(j))) = (j). As time reverses the wave amplification is converted to attenuation (damping in some sense) 
A{<j>,)A{F{4>-:)) = l. 

In unabridged mathematical notations we have the function 

sin [ip — a) cos p + sin [Lp — p) cos a 
defined in the interval Lp S (— 7r/2, 7r/2). For arbitrary values of the parameters a and (3 

F[Fm ^ 0, (83) 
i.e., this function F coincides with its inverse function F~^ . The nonlinear function F has only 2 immovable points 

F{a) = a, F{f3) = /?. (84) 

Defining also 

A{(p) = I [sin {(fi — a) sin /? + sin {(p — /?) sina]^ + [sin (ip — a) cos /? + sin (ip — j3) cosq!]^| / sin^(a — /?), (85) 
we have another curious relation 

A[F{p,)]A{p)^l. (86) 

The so derived amplification coefficient A{(f>i; Sg, S^i) depends on the initial phase. In the next subsection we will 
consider the statistical problem of phase averaging with respect to the initial phase (t>i. 
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C. Phase Averaged Amplification 

For waves generated by turbulence the initial phase is unknown and one can suppose a uniform phase distribution. 
That is why for solving the statistical problem of energy amplification we need to calculate average values with 
respect to the initial phase (p^. That idea is coming from the well-known random phase approximation (RPA) in 
plasma physics. The phase averaging already introduces an element of irreversibility because we already suppose that 
waves are created with random phases. This is the MHD analog of the molecular chaos from the theorem for entropy 
increase in the framework of the kinetic theory if the probability distributions are introduced in the initial conditions 
of the mechanical problem. In the case of accretion flows we also suppose that turbulence is a chaotic phenomenon 
and we have to apply the RPA for investigating the statistical properties. 

The calculation of the integral 

^ ^AA(0O = 2 - s^g = 2 - P (87) 







jives for the initial phase averaged gain 

fTT/2 



J-tt/2 Tl" 



l^l- (88) 



/2 I" sin ((5u - ^g) 

In such a way the SMW amplification coefficient Q is presented by the tunneling coefficient V of the corresponding 
quantum problem. Both coefficients are expressed by the asymptotic phases 5g, iSu in analogy with partial waves 
phase analysis of the quantum mechanical scattering problem in atomic and nuclear physics. The axial symmetry of 

this result Q K"^ + K'^ significantly simplifies the further statistical analysis. 

Let us analyze the physical meaning of the gain coefRcient Q. As 

= KyT = Vp^ky t = LUASgn{ky)t, tJA(k) = |Va • k| = VA\ky\ > 0, (89) 
Vgr = = VASgn (Va • k) = VASgn(fcj^)e.y, (90) 

the asymptotics Eq. (|73p let us conclude that for |t| — + oo we have only magnetosonic waves with dispersion coinciding 
with that of Alfven wavesi^ In the spirit of M. T. Weiss quantum interpretation of the classical Manley-Rowe 
theoreroiii^ one can present the wave energy HujaN by a number of quanta, the number of alfvenons: "The alfvenons 
introduced in this Lcttei— appear to be effective and spectacular converters of electromagnetic energy flux into kinetic 
energy of particles." We use this notion in a slightly different sense, our former terminology was alfvonaiS^ in our case. 
Following this interpretation, the energy gain Q describes the increasing number of quanta 

S (xf^[^ - fi^A(|fc,|)galfvcnon.(t^+00)^ ^g^^ 



^A(|fcal)A^alfvcnons(i ^ -Oo) 

as in a laser system. In this terminology the mechanism of heating of quasars can be phrased, namely due to lasing of 
alfvenons in shear flows of magnetized plasma. Laser or rather maser— effects are typical phenomena in space plasmas. 
The hydrodynamic ovcrreflcction instability^ and burst-like increase of the wave amplitude^^ are phenomena of similar 
kind. More precisely Q is the gain for the z-y-polarized SMWs, the energy of mode conversion in z-polarized AWs 
will be analyzed elsewhere. The notion ampliflcation is correct for standing waves but 2 cos(k • r) = e* '' + e"* '' and 
amplification is simultaneous for opposite wave- vectors. That is why some theorists prefer to use overreflection in 
spite that there is no rigid object reflecting the waves. 



D. Analytical Approximations for Amplification 

For scattering problems by a localized potential at small wave- vectors Q <C 1, when the wavelength is much larger 
than the typical size of the nonzero potential, we can apply the delta-function approximation 

|^t^(o-(rra^-2go<5(e). (92) 

In this well-known quantum mechanical problem^^ the transmission coefRcient is 

"""WTqI^ ^" = "2' ^^=2+"'^^*"^Q^' ^« - cos(J,) ' 
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According to the tradition of the method of potential of zero radius, the coefficient Qo w ^tt is determined by 
the behavior of the phases at small wave-vectors. Only qualitatively this parameter corresponds to the area of the 
potential 

/oo ^ 
^ (1 + ^2)2 ^^^ (94) 

but the renormalizing coefficient Z^cn ~ 2 which we have to introduce differs from unity. According to the general 
relation Eq. ([55)1 for the amplification we have 

c;-l = 2(^i-l)«^, 2Qg = ^« 4.934«5. (95) 



The scattering phases Sg and 5u can be obtained by a fit of the asymptotic wave functions Eq. ((M)) and Eq. 
with the analytical solutions Eq. ([51]) and Eq. ([5^ . In such a way we can calculate the wave- vector dependence of 
the transmission coefficient as it is depicted in Fig. [TJ 

Let us derive the analytical formula for Qq: The solutions 



MO =v^Te, (96) 

MO = v/l+e'arctan(0 (97) 
of effective Schrodinger's equation Eq. ([27)) in long wavelength limit Q ^ 

4^ - (TTa^^ = ° 

with asymptotics at ^ ^ oo 

^3(C»1)«^, (99) 
V'«(C»l)«f^ (100) 

have to be compared with the approximative solutions Eq. ([M)) and Eq. (|65p for ^ ^ 1, when for Q ^ we have 
\^g\ ~ \^u \ ~ f and sinusoids are almost linear 



1 

Q' 
1 



V-gll « e « 7t) « Dg sin(Qe) - DgQ^, (101) 



^„(1 « ^ « g) ~ ^« sin(QO « -DnQe (102) 
The comparison of the first derivatives in this region 

de^g(C » 1) « 1 « DgQ, (103) 

d^Vn(C»l)«|«^nQ (104) 



gives 



2Q 

Then formula for Wronskian Eq. (|7^ determines the phase difference 



^9 « Du - (105) 



TT 

and transmission coefficient 

2 



sin(5g - (5u) « -Q (106) 



V^sin^iSg-S^) = ( ^ ) . (107) 
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-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 

log Q 

FIG. 1: The efficiency of tfie dissipation meclianism that is likely responsible for the formation of stars and other compact 
astrophysical objects: logarithm of the wave energy amplification Q = 2/'D — 1 (cf. with Eq. (|88}) versus the logarithm of the 
dimensionless wave-vector Q. The delta-function approximation Eq. (|95p (upper cur\'(') is asymptotically exact at huge (giant) 
amplifications. In the inset the transmission coefficient D is plotted versus the square of the wave-vector of the auxiliary 
quantum mechanical problem. The delta-function potential approximation (lower curve) Eq. (|93|) is good only for long enough 
wavelengths ai which the amplification is significant 6^1; log 5 = 6 means energy amplification one million times or 

60 dB. 



The Eq. ([55)) then gives for the phase averaged amplification 

2 



G = — - 1 w 7^ + const (108) 



and the comparison with (5-potcntial approach Eq. (|95p gives the analytical result for the strength of the (5-potcntial 
Qo = 7r/2 given in Eq. ([55]) . 

Our problem formally coincides with the quantum probleroi^ of transmission coefficient P cx i? at low energies 



[1 + [ax)'']'' 2m V I" / Uq 

The zero-radius potential Padc approximant Eq. (|93p has an acceptable accuracy for small wave-vectors. Having 
the quantum mechanical transmission coefficient, we can calculate the energy gain coefficient Q — 1. In the next 
section we will incorporate the approximative solution 

in the simplest model of Kolmogorov turbulence. Here the subscript _L means that the amplification depends on the 
projection of the wave-vector perpendicular to the shear velocity and magnetic field. 

The delta function approximation has a visual interpretation in classical mechanics as well, supposing that -0 is the 
displacement of an oscillator and ^ is the time. The approximative equation 

d^^ = -Q'HO + 2QoS{Om (111) 
means that at a time moment ^ = the oscillator is subjected to a forcing impulse with a magnitude 

d{V(+0) - = 2Qotjj{0), or v^.{+0) - w^(-O) = ^'^''^^'> ^ when i^^(O) = 0. (112) 
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If for Q <C Qo we have initial oscillations with amplitude Ai 

TZ-^AiCosQ^ for C<0, (113) 
after the push in ^ = we have oscillations with much increased amplitude 

il; = Ai singe, for ^ > 0, Ai ^ ^A-, > A-,. (114) 

The strong push with an appropriate phase "amplifies" the oscillations; this burst-like increase of the wave amplitude 
was observed in numerical investigations of linearized two-dimensional MHD equationsi ^^i^^ This phenomenon is akin 
to the extremely strong hydrodynamic instabilities due to a velocity jump; its prediction and discovery both in theory 
and experiments are described in Ref . [1^ . In such a way the Alfven's idea of the importance of MHD waves in the 
transfer of momentum is reduced to the very simple mathematics for the jump of the velocity of SMWs Eq. (jll2p . It 
is instructive to rewrite this force in the r-space. 



IV. INCORPORATION OF TURBULENCE AS RANDOM DRIVER OF MHD WAVES 

Being as efficient as has been demonstrated above, the amplification of SMWs, Eq. (jllOp . is possibly the dominant 
physical factor responsible for generating and maintaining the turbulence in accretion fiows. However the theory of 
turbulence is much more complicated than the theory of linearized waves. That is why here we provide an illustration 
how the wave amplification can be incorporated in the turbulence theory. In order to establish common set of notions 
and notations we will recall some basic properties of the homogeneous isotropic Kolmogorov-Obukhov turbulence. 



A. Kolmogorov Turbulence 



V(r) = / e^-Vk-^, Vk = / V(r)c-"'--d3x. (115) 



Let the velocity be presented by the Fourier integral 
The energy per unit mass is 

d^fc f „ d^fc 



where we introduce the spectral density averaged with respect to the turbulence 

= ^(Vj^)turbulonce- (117) 

We can introduce also the energy of vortices which contains Fourier components with wavelength 27r/fc, shorter than 
some fixed length A 



where for isotropic turbulence = £k- This energy evaluates turbulent pulsation with size A. Further on we will 
continue with only order of magnitude evaluations, hence in the following estimations we will drop off factors such as 
47r, i, etc. 

According to the Kolmogorov-Obukhov (KO) scenario in the inertial range the magnitude of the velocity pulsations 
V\ can depend only on the turbulent power dissipated per unit mass e. There is only one combination with the 
appropriate dimension 

energy/mass power 
— ' - (119) 



\/V\ time — length/ velocity mass 
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which yields 

lkX>l JkX>l 1^ ' 



E{k) = J k'S^^dn^ C^e^'^k-^'^ Ck « 1.6, f ^ ^ (120) 

Vx is the amphtude of variation in the velocity pulsation at distance A. £k is the energy density in the k-space per 
unit mass; in the Kolmogorov-Obukhov picture this is a static variable. 

The scaling law V\ ~ (eA)^/'^ Eq. (|119[) is applicable for large enough distances A > Aq, where Aq describes the scale 
where dissipation effects become essential 

71,^1/2215/2 

\^ = {yl/eY'\ y^,^il/p, p = MNp, Np = N,, ^ = 0.4— (121) 



p = NpTp + N^T^, -« 1836, -3 = 3^^ = ^, ^ 'p = e^MpM^^Cp ^ 



e4mi/2/:^ 



< 1, 



eVmi/2/;^ . 1 eW/2£, _ i^k _ 0.24 x 4^ri*/'rp'/' ^ 

J^m ~ EqC — TT^, £> — -; Pm — — r at ^ ^■ 

0.6 X AnTe' 47reo 0.6Te ' ^tic^ eb^Vp 

Here, apropos, we introduced self-explanatory notations for the parameters of fully ionized hydrogen plasma: bare 
plasma viscosity r/, the Coulomb logarithm for protons Cp, the proton mass M, mass density p = MNp, pressure p, 
sound speed Cg, magnetic diffusivity fm, kinematic viscosity i^k, heat conductivity n, specific heat per unit volume 
Cy, temperature conductivity Xi Prandtl number P, magnetic Prandtl number Pm, total viscosity determining Alfvcn 
waves damping i/ = i^k + i^m, Ohmic resistivity g, constant in Coulomb interaction e^, Debye screening length Ad, 
electron temperature Te and proton temperature Tp times the Boltzmann constant, number of electrons and protons 
per unit volume = Np, etc. Those formulas are system invariant: in SI po ^ 4tt x 10^ — I/c^Eq, in Gaussian 
system //g = 47r = 1/eo, while in the Heaviside-Lorentz system po = 1 = l/^o. 

Let us now consider a magnctosonic wave with a time-dependent wave- vector Eq. ^ 

qx{t) — —A{t — to)qy, Qy ~ const, qy ~ const, (122) 

and time-dependent energy density per unit mass in real space 

w{t) = i(VLve + ^> = ^ k + (dQcb)'l , (123) 

Z PPO 4 L J 

where (...) stands for spatial averaging, (cos^(k • r)) = i. Then the energy density in the k-space is 

£^{t) ^ w{t)5[k ^ q{t)]. (124) 
Let us mention that all MHD variables b, v, w in Eq. (|123p . and P in Eq. ([M)) depend on the effective wave functions 



ijj and X (solutions to the effective Sclirodinger equations) through the dimensionless time ^: 



_2Kl f 

KyQj^ 
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(125) 



I MQ{i^O]^^Ai' + x{0. (127) 



r eos[o(^-ni^df' + ^(i±i^lM(i)^ 



(128) 

da(0, (129) 



cos[Q(^-0]T-^de' + ^. (130) 



Q 7-00 ' " "1 



For programming one can use also formulae explicitly expressed by the initial conditions b(to), v(io) 



^ ' Tr^ + Tr(^-^o), if.(r)=/v,(ro)-/s(r-ro), (131) 



2^' r • /^.N. .^.A./.M , Q ^V'(0 ..a _ Q 



= -TTi:^ Nn(gC)/c(0 - cos(gO/.(0] + ^^== - T^xiO = ^^fc. - 17^6., (132) 



Qx/ V-.. .v.. v~.s. .V... yiTe^ X,-^^^^ A,^-^ Ky 

sir 

(i + e')d?V'(e)-ev^(e) 



2K' 

= ^ [sm(Qe)ic(0 - cos(gO/s(0] + X(0, (133) 



where 



Q(l + ^2)3/2 "dQC?^., (134) 

Ay Ay 

- [cos(Qe)ic(0 + sin(ge)/s(e)] + dga. (136) 

m . £ cos(gn ^ d, ^ £ cos(go ^i±«™^d,, (13T) 

. £ sin(gn ^ d, . £ sm(gn ^^iMM^dC, (138) 

x(0 = &.(eo)cos[g(e-eo)]+^'.sin[g(C-eo)], ^0 = (I39) 

dQ5x(e) = -&.(Co) sin[g(e - eo)] + V, cos[g(e - Co)], K^io - -g^ = if.(eo) + g(e - eo), (mo) 

4^{S.a) = Jl + ilh^. d^Hio) = Qjl + ilvM + ^k=bM, (141) 

V 1 + 

V'(0 = [^(Co)dcV^u(Co) - V'u(^o)d^V'(eo)] V-glO + [^g(eo)dcVXeo) - V'(eo)d«^g(^o)] V'u(^). (142) 
In longwavelength limit g <C 1 and ^0 ^ ^00 the numerical integration gives 
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FIG. 2: Radial wave component of the magnetic field bx{£,) as function of the dimensionless time ^. 
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FIG. 3: Time dependenve of the radial velocity Vx{S,)- 



Due to odd integrants Jc,g = = ,/s,u; liniQ^o[Q£'u] = 1.571. 

The formula for bz describes mutual transformation between MHD waves with orthogonal polarization. Only for 
Kz ~ we have exact separation between SMWs with xy-polarization and AWs with z-polarization. 

An illustration for the initial conditions ?/'(— 100) = 1, dj?/;(— 100) = 0, x(~c«) = = djx(— oo) and parameters 
Ky = 0.3, Kz = 0.1, is presented at the figures below. For reliability and check of formulae the figures are doubled 
by the numerical solution of the set Eq. (fT2|) by Runge-Kutta methodi^ Those general formulas give a solution 
to the Cauchy problem. Having in the beginning t = to a distribution of the magnetic field BwaveCr, to) and velocity 
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FIG. 4: Phase portrait of the radial r — motion Vx versus bx. The big loop outlined at ^ = is the driving force of SMW 
amplification. This kick amplifies magnetosonic waves, but the radial amplitude disappears at ^ = 100. The curve starts and 
ends at the point (0,0). 

with V • Vwavc(r,io) = 0; we can calculate the Fourier components 

vu(to) = i / Xz2^^c-'''-dx3, (145) 

bk(io) = / e(r,to) ^_i^.,^^3^ ^^^g^ 

J ^0 

Co,k^-^i=, (147) 

and initial dimensionless time ^o,k- li ~ then sgn(A;j^)^o. k = '''o.k = —k^/ky. Then we have to determine the 
coefficients C in the general solutions for ip Eq. ([^5]) and x Eq. (HH) using the initial values at to 

bk(e) = Cgbg + CuK + Cgbg + Cuba, (148) 
Vk(6 = CgVg + ab„ + Cgvg + Cuva, (149) 
k-bk = = k-vk. (150) 

In this set we can use only x- and z-components, and so we obtain 4 equation for the constants Cg, Cu, Cg, and 
Cu. The functions bg(^) and Vg(^) are defined via substituting Vg in Eqs. (fT25l) - (fT30)l or Eqs. p^ - pl2l) . and 
analogously ipu, Xg, and Xu- Then at each moment t we can calculate all variables in the k-space 

aw =eo,k + (t-to)A^==^, (151) 

bk(t) = bk(ek(i)), (152) 
Mt) = Vk(ek(i)), (153) 
fcr™(^) ^k^-{t- to)Aky. (154) 
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FIG. 7: Phase plot Vy versus by of axial — SMW oscillations parallel to the constant magnetic field Boip. The wave 
amplification is the ratio of the cycle areas at ^ ^ oo and ^ — oo. The fast transition between those two orbits describes the 
lazing of the alfvenons. 



Finally, we can return back to the real r-space 

Vwavo(r,t) = -i^A J^k{t)c 
Bwavo(r,t) = -Bo y"bk(t)c 



-i[]i-r—{t — tQ)Akyx] 



Re(e" 



cos if, 



i[k-r— {t — to) Akyx] 

Re(-ie-''^) : 



(27r)3' 

dfc3 



(27r)3^ 
simp. 



(155) 



This evolution of MHD variables is the main detail of the theory of MHD turbulence in a shear flow. 

Consider now an imaginary fluid filling the phase space k and w(i) = £'k(i) from Eq. p24p being the energy density 
carried by a droplet of that fluid. As a wave mode initially with wave- vector k evolves according to Eq. (|122p . the 
infinitesimal phase-fluid droplet associated with that mode moves in the k-space. Wave amplification means that the 
energy density of the droplets increases by a factor of G 



G = 



w{t +oo) 
w{t — oo) 



Indeed for x — 0, kz = 0, and ^ ^ oo 



Q 



and 



w{t^^) = -VlDl 



(156) 



(157) 



(158) 



For big enough time arguments |^| 3> 1 and purely two-dimensional waves with fc^ = the motion of the fluid 
asymptotically corresponds to a SMW with dispersion coinciding with the AW one 



ipiO ^ A cos(w3„„t -f 0f). 



(159) 
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FIG. 8: Time dependence of the magnetic field oscillations bz{£,) perpendicular to the disk plane. 
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FIG. 9: Out of plane velocity oscillations parallel to the disk axis z as function of the dimensionless time Vz{^). 

The Poynting vector, i.e., the energy flux in r-space is Vaw. 

The velocity of the droplet in the k-spacc according to Eq. (|122p determines the field of the shear flow in the k-space 

\J = dtq{t) ^ -Aqye.,, Vt''^' = -Akye,. (160) 

Looking at a droplet we actually derive the shear flow velocity field in k-space, U^'^'''*''. 

According to the Kolniogorov-Obukhov cascade of energy we have a constant energy flux through each spherical 
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FIG. 10: Phase portrait of the AW oscillations transversal to the disk plane; Vz versus b^. The AW polarization is almost 
perpendicular the constant component of the magnetic field and the wavevector for ^ = 100. The wave amplitude at ^ = —100 
is negligible; we have mode conversion of SMW to AW at ^ = 0. 
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FIG. 11: Dimensionles wave energy density w = i(v^+b^) versus dimensionless time ^. In the inset the ordinate is logarithmic. 
The ratio ii;(100)/?i)(— 100) describes the energy amplification by the shear flow which is the heating mechanism of accretion 
disks. The analytical expression of the energy gain Eqs. (|13H4142]) is given by the confluent Heun functions Eqs. (|31II32|I . The 
approximative central symmetry of the blue curve corresponds to 10% accuracy of the product ■w{t)w{—t) ~ const for this 
numerical example. 
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surface with surface element df in k-space 

e-^fj^OuJ^Odf = 4KOc/feK°4^fc2, (161) 



which gives 



Uj^o^eV3fcV3e,, efe = ^, (162) 



i.e., the velocity in k-space has dimension 1/ (time x length). Here we used an important for our further work notion 
of the energy flux in the k-space 

S = fkUk (163) 

which is equal to energy density times velocity in the k-spacc. This notion is analogous to the Poynting vector being, 
however, defined in the k-space. In the Kolmogorov-Obukhov scenario we have 

A ■ S^o = sS{k). (164) 

In order to approximate the turbulence as an initial source of MHD waves we have to merge the turbulence with 
the wave spectral densities and velocities. The simplest possible scenario is given in the next subsection. 



B. Derivation of Shakura— Sunyaev Phenomenology in 
the Framework of Kolmogorov Turbulence 

How vortices create waves is a complicated problem far beyond the scope of the present study. Here we will give 
only a model illustration merging the spectral density of vortices E^^'^^ from Kolmogorov turbulence with spectral 
density of magnctosonic waves 



on the plane in momentum space 



sgn{ky)A-\ kl + kl<A-~^, (166) 



where we qualitatively suppose that vortices are converted into waves. Sign function corresponds to the direction of 
the shear flow in the k-space, Eq. (|160p . For ky we consider that turbulent vortices have a given spectral density at 
kx > A^^ which is converted to MHD wave energy at the plane k^ = A^^, and further on this wave energy evolves 
according to our solution. In other words, the plane kx > A~^ is the boundary between the vortex region and the 
beginning of the amplification in the wave region where vortices have negligible influence. In our qualitative picture 
we suppose that vortices create spectral density which further on evolves as wave spectral density with negligible 
influence. 

The amplification Eq. (jllOp is essential Q ^ I only within a cylinder 

Giky, fc,) - 1 ^ = kl + kl < A-2 (167) 

with radius A^^. This result with remains imchangcd in amplitude if we include the Jc,u and Jg^g terms. 

The amplification occurs in the region — A^^ < k^ < A^^, that is to say from the cylinder we cut a tube with length 
2A~"'^. In other words, we have a domain with a shape of a tube in momentum space 

V^{kl + kl<A-', |fc,|<A-i}. (168) 

In order to calculate the total power of waves H (per unit mass) analogously to Eq. (|16ip we will integrate the energy 
flux on the surface of the tube 

J^'^'^^-'df. (169) 
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As the shear in the physical flow resuhs in a drift of the wave modes along the axis of the tube, we have to take 
into account only the circular surfaces 

e - J |C/x| [g{ky,k,) - 1] S^dkydk,. (170) 

The multiplier [Q — 1) takes into account the difference between flowing out and flowing in energy fluxes. 
We can use polar coordinates 

kz—q cos 9, ky — q sin 9, Ux — Aqsin9. (171) 

Averaging over the angle 9 

2 3.9 2 

(Ux) = -Aq - Aq, / sin^— = -, 
TT n n 

and substituting it in Eq. (|169|1 . using dkydkz = d{7rq'^), leads to the simple integral 



(172) 



^SMl-^-AVl (173) 

Then for the volume density of the amplifled waves we have 

Q = pH^ pAVi ^ pieVA)^/^A^/\ (174) 

As all waves are finally dissipated, Q is actually the volume density of plasma heating. 
For evanescent Kolmogorov turbulence power 

pe « ApVi = AB^/nn (175) 

the heating power H has a critical behavior 

Ti / AV^ \ 

d,n - Gturb = 7 ( ) > 1' ^ > e for e -> (176) 

which demonstrates that disks can ignite as a star even for very weak turbulence and magnetic field. The ratio of 
wave power and Kolmogorov vortex power Gturb can be considered as an amplification coefficient for the turbulence. 
This energy gain shows how efficient is the transformation of shear flow energy into waves or in a broader framework 
the transformation of gravitational energy into heat of accretion disks. 

For hydrogen plasma pc^/p = 5/3 '--^ 1. Now we can evaluate the shear stress (as given by the ratio of the volume 
density of heating power and the shear frequency) 

" = :4-^n — ) ^^^^ 



via an effective viscosity 



and the dimensionless Shakura-Sunyaev friction coefficient 

(179) 

Including of the energy of z-polarized AWs does not modify this result. Here we wish to emphasize that in our 
evaluation of the energy gain, we were concentrated on the wave amplification of the energy of two dimensional 
motion in the x-y plane. Taking into account the energy in z-direction shows that the heating is even higher, which 
is of course in the favor of the concepts. 

For an approximately Keplerian disk rotation the shear rate is half of the frequency of the orbital Keplerian angular 
velocity A = — ^WKcpior- In this case for time A~^ the disk rotates per 2 radians. For Earth's rotation along the 
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Sun this time is of the order of one season. In such a way the length parameter of our problem K = Va/A can be 
evaluated as one Alfven season. Then Va from the final result for the Shakura-Sunyaev parameter can be qualitatively 
considered as a pulsation of the turbulent velocity for two disk particles at distance equal to one Alfven season A. 
Our theory is formally applicable for Va ^ Cg but the boundary of its applicability (when compressibility effects 
stop the SMWs amplification) allows us to understand that strong disk's turbulence can lead to Shakura-Sunyaev 
upper limit a ~ 1. Thus the following cascade of events emerges as a likely scenario for the intense heating in 
accretion flows: the heating of the bulk of the disk creates convection. For strong heating the convection is turbulent. 
Turbulence generates magnctohydrodynamic waves. Waves arc amplified by the shear flow - this is the transformation 
of gravitational energy of orbiting plasma into waves. Waves finally arc absorbed by the viscosity which creates the 
heating. The heat is emitted through the surface of the disk. This process of formation of stars and other compact 
astrophysical objects from nebulas works continuously - we have a self-consistent theory for self-sustained turbulence 
of the magnetized accretion disks. 

The weak point of this scenario is the supposed convective turbulence which in presence of magnetic fields is 
unlikely to be of Kolmogorov type. We consider as much more plausible scenario the appearance of a self-sustained 
magneto-hydrodynamical turbulence considered in the next subsection. 



C. Kraichnan Turbulence as a more Plausible Scenario for Accretion Disks 



Magnetic field qualitatively changes the behavior of the fluid. Wc have no waves generated by vortices - the 
turbulence in magnetic field is related to MHD waves. Analogously to the Kolmogorov law Eq. (|119p . for the Kraichnan 
turbulence the power of energy cascade in the dissipation-free regime is given by the wave-wave interaction 

„ _ iVx? _ (velocity)^ _ power 
AVa length mass 

This power is proportional to the intensity of the two interacting waves and this nonlinear effect for incompressible 
fluid is due to the convective term V • VV of the substantial acceleration DfV = dtV + V • VV of the momentum 
equation Eq. 

The theory of generation of SMWs invokes parallels with the nonlinear optical phenomena in lasers. The velocity 
oscillations of two amplified MHD waves V^'^' and V(^) create an external driving force of the new wave with velocity 
field V. In the linearized Eq. ([8]) we have to insert a small nonlinear correction 

pdtV = -Vp + ^ ^ X B + pf , f = i Vv*'*' • VVt*"'. (181) 
^J■o 2 ^ 

a, f> 

Here, in the inhomogencous term f we have to perform summation over all other MHD waves. This external for the 
wave force (per unit mass) acts as an external noise and its statistical properties are determined by the force-force 
correlator 

f(ri,ti;r2,i2) = (f (ri, ti) f (ra, ia)), (182) 

where the averaging is over the waves phases. A scenario of such type (a Langevin MHD) was described in Ref. [l^; 
this approach is similar in the spirit to the forced burgers turbulencei^ In the framework of that scenario the strongly 
amplified \D{\ ^ 1 MHD waves with asymptotics Eq. (|73p 

ij^Df 9{^) cos(Q^ + 0f ) exp {-i^'K^rye) (183) 

generate new waves and after a statistical averaging we have a self-consistent theory for magnetic turbulence in a 
shear flow. So MHD waves ignite the chain reaction of quasar self-heating. The last exponential term describes the 
wave damping when a small viscosity is taken into account. Damping is significant only for t ^ oo when \kx\ 2> \ky\ 
and the wave-vector is almost parallel to the magnetic field. In this geometry, the damping rate of the wave density 
of AWs and SMWs is proportional to the square of the frequencyii 

w{t) = w{0) exp (^-^lyt^ = w{0) exp {-vk^t) . (184) 

For the time-dependent wave- vector fc^(t) w (kyAt)'^ in the argument of the exponent we have to make the replacement 

i^kh ^ u f e{t')<lt' = \i^'KIt\ = ^ = 1, 7^ ^ 5mri ^ (185) 
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where v' is the dimensionless viscosity and Ti. = Vp^/A is the "Reynolds number of the magnetorotational instabihty 
(MRI)^" and 5mri is the Lundquist number of MRI. After long enough time when 

|fc.(t.)| = 1/Aa, Aa=T^, (186) 

Ka 

MHD waves are completely dissipated. 

The details of self-consistent MHD turbulence will be given elsewhere, but again the wave amplification operates 
as a turbulence amplifier. For MHD turbulence one can expect 

o-_R<p = a,^{v')pB, Pb = \p^A- (187) 

The evaluation of the magnetic friction coefficient am as a function of the dimensionless viscosity is a new problem 
addressed to the theoretical astrophysics. 



V. DISCUSSIONS, CONCLUSIONS, AND PERSPECTIVES 



We propose that the long-missing element of the dissipation mechanism of accretion disks is now identified - it 
is the amplification of long-wavelength SMWs. This is the indispensable ingredient of the stars generating engine. 
Without it the Universe would possibly be a structureless gas - deserted and uninhabited. The two-dimensional 
disks and the redistribution of angular momentum in them are the means of creating a one-dimensional compact 
astrophysical object from gases and dust. The planetary system is the result of this star-producing sequence of events 
when the accretion stops and the disk is frozen. It is remarkable that a simple equation of the Schrodinger type, 
Eq. ((27)) . is at the core of the friction mechanism which cheated the diversity of the Universe. This means that the 
Schrodinger equation can describe one more phenomenon while science is on its way to explaining the frogs and the 
musical composers. "The next great era of awakening of human intellect may well produce a method of understanding 
the qualitative content of equations. Today we cannot. Today we cannot see that the water flow equations contain 
such things as the barber pole structure of turbulence that one sees between rotating cylinders. Today we cannot see 
whether Schrodinger's equation contains frogs, musical composers, or morality ~ or whether it does not.—" 

Notably, here we have observed only an amplification in a weak magnetic field but not instability, where the 
amplitude of waves increases infinitely with time. Also, the rotation is found not relevant for this phenomenon. 
Therefore, regardless of similarity in spirit, our work is completely different from the research focused on rotational 
instabilities in strong magnetic fields. Our results point to an amplification but not to an instability. For sufficiently 
long wavelengths the amplification can be enormous but never infinite. Amplification means exponential increasing 
of the amplitude in the framework of linear theory. 

After investigating the local dissipation and shear tension by the self-consistent statistical MHD method the cor- 
responding numerical value of the a parameter can be incorporated in global models for accretion disks dynamics. 
The global models for accretion disks include also the problem of disk dynamoi^ In order to create a self-consistent 
magnetic field i?o,v l^e accretion disks operate as radial inflow generators from plasma physics;^ see also Ref. HO. 
In this broad program the local and detailed investigation of magnetosonic waves propagation in a homogeneous 
shear flow and magnetic field is only the first step in our understanding of the properties of space plasmas, i.e., the 
low-gradient approximation is the indispensable step in our understanding of accretion power in the Universci^ 

The further development of the theory of MHD waves turbulence will give additional important details but even 
with what we know now, we get an insight into the workings of the star-creating engine. A very powerful, hence 
possibly dominant, mechanism of energy transformation in shear magnetohydrodynamic flows is identified, which 
is the amplification of Alfven waves. They are also responsible for the heating of the solar corona— and represent 
a quite common mechanism for heating of space plasmas in general. Here we wish to recall that shear flows are 
important for formation of toroidal magnetic flcld from poloidal one. The differential rotation^S of the Sun is another 
shear flow which can lead to ampliflcation of MHD waves. The convection excites MHD waves in almost toroidal 
magnetic fluxtubcs. When buoyant magnetic fluxtubes reach photosphcric surface the MHD waves can be significantly 
amplified. 

Recently Alfven waves of sufficient strength have been unambiguously observed in the chromosphere by Solar Optical 
Telescope onboard the Japanese Hinode satellite. Such Alfven waves are energetic enough to accelerate the solar wind 
and to heat the quiet coronai^ In order to reach quantitative agreement of the theory it is necessary to merge the 
spectral density oc Daw/uj" of these AW with the speed of the solar wind Wwind and coronal temperature Tp. The 
mission of the theoretical models is to give the simple relations between those experimentally accessible parameters 
/ (Tp, Wwind, ^Aw) = Oi cf. Rcf. (ssj . Heating of tokamak plasmas by magnetosonic waves is also a widely discussed 
issuei^i^ According our scenario the self-heated accretion discs are just working tokamaks. 



26 



The idea of accretion disks can be traced in the development of contemporary sciencci^ While for Descartes (1644), 
Kant, and Laplace that idea was a creative mythology;^ the Hubble Space Telescope has now delivered observational 
evidence for disks of ionized gases around massive black holes.— Observations of protoplanetary disks are discussed 
at the web-page^ and its linksi^i^i^ 

It was only several decades after a horrific punishment for the expression of now mainstream views took place at the 
market place Campo dei Fiori^ in Rome, when Galileo observed spots on the Sun and detected the Sun's rotation. 
Why the Sun is rotating so slowly (i.e., the problem of angular momentum dissipation) eventually becomes a focus 
issue in cosmogony 4^ The development of plasma physics determined that the molecular viscosity of plasma is far too 
small to account for the observed heating and that new ideas should be tested. The amplification of SMWs, an idea 
with the potential of explaining the enormous luminosity of quasars, was explored in the present work. 

Pursuing further applications and proof of concepts, one could notice that the magnetosonic waves of magnetic 
turbulence can emit radio waves through the large surface of the accretion disk. Correlation between the radio waves 
and optical emission from quasars will be the crucial test for the present theory of dissipation in accretion disks. 
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